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SUMMARY 

A methodology has been developed for the computational simulation of struc- 
tural fracture In fiber composites. This methodology °* 5 | J p 

Drocedures for mixed mode fracture in generic components and of an ’ntegratea 
comDuter code CODSTRAN (Composite Durability Structural Analysis). The ge e 
"es of composite structural fracture include:. (1) single and combined mod 
fracture in beams, (2) laminate free-edge delamination fracture i and <3> lami- 
nate center flaw progressive fracture. Structural fracture is ide^ti ^ 

?apid chan es in one or all of the following: (1) displacements, 2) frequen- 

cies (3) the buckling loads, or (4) the strain energy release rate. These 
rapid changes are herein assumed to denote imminent structural fra cture. ase 
on P these rapid changes, parameters are identified which can be used as guide- 
lines for (1) structural fracture, (2) inspection intervals, and (3) retirement 

for cause. 


INTRODUCTION 


It is generally accepted that flawed structures fail when the flaws grow 
or coalesce to a critical dimension such that (1) the structure . cannot sa y 
perform as designed and qualified or (2) catastrophic fracture is imminent. 

This is true for structures either made from traditional homogeneous materia 
or fiber composites. One difference between fiber compos i tes an rad t ona 
materials is that composites have multiple fracture modes that in tiate local 
m»s Spared to only a few for traditional materials Any predic ive app oaoh 
to simulate structural fracture in fiber composites needs to Borina 11 y quant If y . 
(1) these multiple fracture modes, (2) the types of flaws they imt , 

(3) the coalescing and propagation of these flaws to critical dimensions o 

imminent structure fracture. 

An ongoing research activity at NASA Lewis Research Center i-s directed 
toward the development of a methodology for the computational simulation of 
structural fracture in fiber composites. A part of this ^^odoiogy cons s 
of step-by-step procedures to simulate individual and mixed mode. fracture 
variety of generic composite components. Another part has been incorpo ate 
into an integrated computer code identified as CODSTRAN for Composi e 
Structural Analysis. The objective of the proposed paper 
fundamental aspects of this methodology and to illustrate 
variety of generic composite structures. 


Durabi 1 1 ty 
is to describe the 
its appl i cation to a 
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The generic types of composite structural fracture include- m „ni a 

Ir^rapld v' n (2! e th: i" ° f *"* f ° ?l 0 *"« = <1> »• J'ipU cZ ™ " 
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In the present study, computational simulation is defined in a ^nerifir 
way. Also general remarks are included with reso t o 1 1 i f 
the procedure to large structures and/or struc?S?Sl sl?ul TD e 

£S,£ HSr - Waw 


COMPUTATIONAL SIMULATION-DEFINITION 

.vokSM a^X^n^c^^ has 

from fundamental concepts, through the use of computers to obU?n the 
desired structural response with acceptable engineering accuracy " 

7 1 ml ted to? ° f the ab0Ve defin ' t,on - computational simulation Is not 

(1) Applied Mathematics 

(2) Approximate Analysis 

(3) Numerical Analysis 

(4) Numerical Testbed 

(5) Solution Algorithm 

However, It may include some or all of the above, plus many others. 


END-NOTCH BEAM FRACTURE 

In this section we describe the application of the structural frarfnro 

i:rlu™^e??„ r %^> aC f^’ f0r b 7" e “'^ P ^°c- dU ^'pe f ^o^?r?a5?e ta ^e? , 
gressive 4 ? e H^rsJbst^lu^r^ ^'^^ei^S'tato^n^e^ 

<«/»•«. t™ 0 ?at?o b <FW?' fr ° m AS - 9raph,te f'lf/epoxy matrix 
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Mode I 


A schematic depicting the strain-energy release rate (SERR) in a double 
cantilever is shown in figure 2. Two curves are shown in this figure. < > the 
solid line for global SERR is determined using the displacement at the point 
where the load is applied, and (2) the local SERR is determined by the crack 
closure technique. For this case, structural fracture is imminent when the 
crack length progresses to about 1.15 in. This length is determined by the 
intercept of the tangents near the rapidly increasing portion of the SERR curve 
(fig. 1). The important observation is that the general concept is applicable 
to even this relatively simple structure. 


Mode II 

A schematic of the end-notches inducing Mode II and combined Modes land 
II Is shown in figure 3. The finite element model and the local progressive 
substructuring are shown in figures 4 and 5, respectively. The details are 
described In Reference 1. Typical results obtained for Mode II are s own in 
figure 6 where the measured range (ref. 2) is also shown. JL 

imminent structure fracture is about 1.18 in. based on the global SERR curve. 

Three points are worth noting: 

(1) Mode II fracture exhibited some stable growth. 

(2) The local approach predicts optimistic results relative to "critical 
crack length but conservative results relative to the critical SERR magni tude. 
The reason is that the global approach incorporates the overall readjustment of 
the structure while the local does not. 

(3) The SERR (G) is within the measured range (ref. 2). 

The Important observation is that the computational simulation procedure cap- 
tures the whole history of the process which leads to structural fracture 
induced by Mode II fracture. 


Combined Modes I and II 

Typical results for combined modes I and II fracture in end notch beams 
are shown in figure 7. Structural fracture is imminent when the crack length 
progresses to about 1.14 in. (global SERR curve). The local curve, again, pre- 
dicts optimistic magnitudes for critical crack length but conservative magni 
tudes for the critical SERR. 

The decomposition to Modes I and II is also shown. These curves were 
determined using the local closure technique (ref. 1). Interestingly struc- 
tural fracture for this condition is driven by Mode I. Also, the curves for 
Modes I and II resemble their respective independent parts in figures 2 ana t>. 
Furthermore, superposition of Modes I and II appears to apply to the local 
mixed mode curve but not to the global curve. The important observation is 
that structural fracture in beams subjected to combined (mixed) modes I and li 
fracture can be computationally simulated and the respective "critical param- 
eters quantified from the general procedure in figure 1. 
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Mixed Modes I, II, and III 

b ^n ‘ e * 1 a rglme ^ S o?%he°r^g I P ^ll?ed ’ ?„ ^ 


ength for various _ unsymmetri c laminate angles are shown in figure 10 
important observations to be made is that the Mode III contribution to 
tural fracture is relatively insignificant compared to the other two 

:p.sL,s: sisysw.ssssr 

cfrr?ic^r re not easi,y or ' - * 


The 
struc- 
A1 so, 
any 


Collectively the four different types of structural fracture described 
JespecTto: structura1 fracture can be computationally simulated with d 


(1) Its history from initiation 
propagation; 


to the onset of rapid (unstable) 


(2) Its critical parameters (for imminent rapid propagation); 

(3) Its predominant fracture modes that drive it 


The very important point to be noted is that all of the above were obtain 

of that - 


LAMINATE FREE-EDGE DELAMINATION FRACTURE 

In this section we summarize the application of the structural fracture to 
components subjected to free-edge delami nations under 1 n-p line load i ng So J 
cif.cally, we describe the following types of delaminations cente^norke? 
and internal. Detailed descriptions on the simulation pro edu"l ^eps are 

Tienvtscrlbed 4 \"M h A I S0 ‘ft Cff6CtS of on Jck ,Sg will 

free edn s hi ndur. J r “ h ?!S t, S ,’"^‘<"9 the origin of the interlaminar 

Utl ? ™> r fs SU foneI h0 “ n f, 9 Ure r ' 2 " here the 3'obal strIVn energy rifease 

1‘ ' eVe,S ° ff by th * ,0 P erc ent°de lamination and 
settles into a stable delamination growth up to 70 percent. 


Center Delamination 

(r P f A fi? h ^H 9 I aPh h° f f C_S S an ° f a ,ab ° r atory specimen with edge delamination 
(ref. 6) and a schematic of its corresponding one obtained by computational 
simulation are shown in figure 13 (ref. 5). The corresponding SERR foJ three 
different composite systems is shown in figure 14. The behavior of a?1 th?ee 

t0 that ln f,9ura Htl'Iof^d 


4 



The observations to be made are: (1) free-edge de aim nation in the pres- 

ence of In-plane loads is benign; (2) different materials «nte readily evalu- 
ated with respect to their structural fracture res ! st !" ce ;„ a "? '^st^ess * 
area may delaminate simultaneously with no change in the in-plane stress 

states . 


Pocket Delaminations 


follows: first, transpiy 


These delaminations are believed to occur as . 

splits develop simultaneously at several places along edge ; 

local delaminations develop in the vicinity of these splits due 9 

stresses. These local delaminations are called, herein, pocket delamina 
tions. Upon subsequent loading, these pocket del ami nations merge 
delaminations . 


to center 


Computational simu 
in figure 15 for three 
delaminations to center 
curve followed by a rap 
for center delamination 
that computational simu 
nation. The "jump" and 
tegic experiments to ve 
interpret "stick-slip" 
bonded joints. 


lation results for these types of del ami nations are shown 
different composite systems. The merging of P° cR et 
delamination is identified by the "jump" in. the SERR 
id decrease and subsequent stable growth similar to that 
in the previous section. The important observation is 
lation captures the progression of this complex delami- 
subsequent "rapid decrease" can be used to guide sca- 
rify this sequence of events. They can also be used .o 
type of progressive fracture prevalent in adhesively 


Interior/Center Delamination 

These types of delaminations may result from lack of bonding. during fabri- 
cation- or thev may occur as a result of inadvertent damage. Their growth to 
imminent structural fracture can be computational ^ simulated a * 4 * ° tensn^ 
of delaminations. Results from these types of simulations for n-plane tensiie 
stress are shown in figure 16 (ref. 5). As can be seen, their influence is 
negligible until large areas (greater than 50 percent) hav ® 6_ 

upon the delaminations have reached the free edges. From here on, t 
like a center-edge delamination. 


The very important point to be noted 
benign in-plane tensile stress fields. 


1 s 


that internal delaminations are 


Delamination Effects on Structural Response 

The delamination effects on structural response can be computationally 
simulated as descried in reference 5. Typical results for the effec 
center-edge delamination on buckling load are snown in figure 17 for three dir 


ferent composite systems. The 
and is relatively small (about 
observation to be made is that 
buckling loads. Although resu 
nations have even less effect. 


decrease in buckling load appears to be linear 
20 percent for 7-percent delaminated area). The 
local delaminations have negligible effect on 
ts are now shown here, i nternal /embedded delami- 
The same observations hold for delamination 


effects on vibration frequencies 
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CENTER-FLAW PROGRESSIVE FRACTURE 


Flaws in plate/shel 1 
nent structural fracture, 
ture at Lewis during the 1 
the CODSTRAN computer code 
Analysis) has been specifi 
progressive fracture in fi 
ing CODSTRAN are: (1) the 

al 1 other modules ; (2) the 
ite mechanics module; and 


type structures propagate to critical sizes for immi- 
Research activities on progressive composite frac- 
ast 10 years have culminated in the development of 
(ref. 7). CODSTRAN (Composite Durability Structural 
cally developed for the computational simulation of 
ber composites (refs. 7 to 9). The modules compris- 
executive module, containing communication links to 
I/O module; (3) the analysis module; (4) The compos- 
(5) the fracture mechanics module. 


A schematic showing the logic of the code 
ical results is shown in figure 18 (ref. 8) 
flaw propagation up to structural fracture. S 
for different flaw types is shown in figure 18 
ditions in figure 19. The structural degradat 
quencies and buckling load is shown in figure 
results from figures 18 to 20 demonstrate that 
ites can be computationally simulated: (1) to 

loading condition, and (3) for any structural 


with an analysis model and typ- 
The code tracks the damage growth 
imulation of progressive fracture 
, and for different loading con- 
ion in terms of vibration fre- 
20 (ref. 9). The collective 
progressive fracture in compos- 
any level of detail, (2) for any 
response . 


GENERAL REMARKS 


Based on the continuing research effort 
are appropriate. These remarks include: (1) 
gained, and (3) future research needs. 


described, some general remarks 
generalization, (2) experience 


General ization 


al 


The computational simulation procedures described previously can be gener- 
ized for structural fracture in composites as follows: 


(1) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


Develop global structural /stress analysis model 
Apply specified loading conditions 
Identify hot spots for these loading conditions 
Introduce flaws in the hot spots 

With specified loads on the global model propagate flaws 
Monitor structural performance degradation (displacement, buck- 
loads, frequencies) versus flaw propagation 

(7) Identify flaw size for unacceptable performance degradation level 

(8) Set qualification, inspection intervals, and reti rement-for-cause cri- 
eria based on the simulated history of structural performance degrada- 
tion versus flaw propagation. 


i ng 


Experience Gained 

The experience gained during this rather extensive research effort lead 
f essence of a suitable environment in order to successfully conduct 
se y defined ong-term research. This research environment includes the 
following essential features: 
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(1) Continuity in research activity 

(2) Collection of participants' with composite knowledge in: 

(i) Structural mechanics principles 
(i i) Finite element analysis 
(i i i ) Composite mechanics 
(iv) Fracture mechanics concepts 

(3) Participants willing to question traditional approaches 

(4) Participants willing to adopt/invent new approaches 

(5) Continuing and unwavering management support 

(6) Availability of computer facilities and supporting personnel 

These features can serve as guidelines to conduct long-term research in compu- 
tational simulation in general. 


Future Research Needs 

Assuming the state-of-the-art on composite structural fracture is as 
described herein, the authors consider the following as near-future research 
needs : 

(1) Incorporation of uncertainties in the simulation 

(2) Development of specialty: (1) finite elements, (2) boundary elements, 

and (3) functions, all of which include local details and capable. or 
capturing steep gradients and thereby provide computational efficiency 

(3) Efficient self-adaptive global /mul ti local scale methods _ 

(4) Time dependence and multiline scaling (local events simulated in 
different time scales compared to global) 

(5) Formulation/programming for parallel processors 

(6) Adaptation to smart structures and health monitoring systems 


CONCLUSIONS 

Based on the results of the research on computational simulation for 
structural fracture in composite the following conclusions are made: 

1. Computational simulation of structural fracture in composite structures 
is ready for extension to wide use applications. 

2. Reluctance to adopt it is natural because of: (i> unfamiliarity and 

(ii) attachment to traditional approaches. 

3. Reluctance can be overcome by: (i) education, (ii) using it. first in 

parallel with traditional approaches, and (iii) make it a specification 
requirement by the procuring agencies in new designs. 
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Behavior of 


• Df.TERMINE REQUISITE COMPOSITE PROPERTIES AT DESIRED 
USING COMPOSITE MICROMECIIANICS AND LAMINATE THEORY 


CONDITIONS 


RUN J-D FINITE ELEMENT ANALYSIS FOR AN ARBITRARY LOAD 


, Lt UAD UNTIL ONE OF LOCAL STRESSES AT ELEMENT(S) NEXT TO 
CRACK TIP REACHES MS CORRESPONDING STRENGTH 


• WITH SCALED LOAD EXTEND CRACK AND PLOT STRAIN ENERC 
VS CRACK LENGTH 


' SELECT CRITICAL "G w 
and CRITICAL V' 

FROM THE INTERSECTION 
OF THE TANGENTS 



RELEASE RATE 

TANGENTS 

INTERSECTION 


• METHOD HAS VERSAMI IT Y/GENERA1 ITY EXTENDED CRACK 

LENGTH, a 

TIGURE 1 . - GENERAL PROCEDURE FOR PREDICTING COMPOSITE STRUCTURES 
FRACTURE TOUGHNESS. 
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FIGURE 2. - DOUBLE -CANT I LEVER ENERGY RELEASED RATE 
COMPARISONS. 



MIXED-MODE-FLEXURE — MIXED MODE (I &II> 

FIGURE 3. - SCHEMATIC OF FLEXURAL TEST FOR INTERLAMINAR FRACTURE MODE TOUGHNESS. 
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ENERGY RELEASE RATE. G, pSl-lN. 



FIGURE 6. - END-NOTCH-FLEXURE ENERGY RELEASE RATE COMPARISONS. 



(SINGLE POINT CONSTRAINED (LEWIS) METHOD). 
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IE MAT 1C OF F.E. MODEL THROUGH-TIIE-THICKNESS DETAILS 
ESSIVE SUBSTRUCTURING). 


2 




OWNING MODE (I) 


SHEAR MODE (II) 


ej 

uu 


oc 




CRACK OPENING LENGTH "A". 


TEARING MODE (III) 


MIXED MODE 


FIGURE TO - firms 01 CRACK LENGIH AND PLY OR I EN I AT ION ON STRAIN ENERGY 
2 ' 



LAMINATE 


ply INTERLAMINAR 

(TREE-EDGE ) 


- LAMINATION DECOMPOSITION FOR FREE-EDGE INTERLAMINAR STRESSES. 


FIGURE IT. 





FIGURE 12. - STRAIN ENERGY RELEASE RATE FOR A 'TYPICAL- LAB 
SPECIMEN (AS/E l±30/90I s 0.6 FVR COMPOSITE). 



87 ksi 
RADIOGRAPH 
[±<l5/0 2 /90 2 ] 


COMPUTATIONALLY 

SIMULATED 

[±30/90 2 3 s 


FIGURE 13. - FREE-EDGE PROGRESSIVE DELAMINATION 
CAN BE COMPUTATIONALLY SIMULATED. 
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LB- IN. /IN. 



DELAMINATION) . 




IQ 20 30 40 50 60 70 


X AREA DELAMINATED 

FIGURE 15. - STRAIN ENERGY RELEASE RATE (6-PLY CENTER/POCKET 
DELAMINATION). 


15 






FRACTURE LOAD 

FIGURE 18. - COMPUTATIONAL SIMULATION OF PROGRESSIVE FRACTURE IN FIBER COMPOSITES. 
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80X MAX. LOAD 


99X MAX. LOAD 


AXIAL TENSION 



87. SX MAX. LOAD 100X MAX. LOAD 


IN-PLANE SHEAR/BENDING 



89X MAX. LOAD 100X max. LOAD 

COMBINED AXIAL TENSION WITH JN-PLANE SHEAR 
FIGURE 19. - PROGRESSIVE FRACTURE FOR INDIVIDUAL AND COMBINED LOADING. 



FIGURE 20. - THE EFFECT OF STRUCTURAL DEGRADATION ON NATURAL 
FREQUENCIES AND BUCKLING LOAD. COMPOSITE: T-3O0/EP0XY. 
(SUPPORT COND. 1 - EDGES PARALLEL TO LOADING FREE) 
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